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We calculate for a binary mixture of Lennard-Jones particles the time 
dependence of the solution of the mode-coupling equations in which the full 
wave vector dependence is taken into account. In addition we also take into 
account the short time dynamics, which we model with a simple memory ker- 
nel. We find that the so obtained solution agrees very well with the time and 
wave vector dependence of the coherent and incoherent intermediate scatter- 
ing functions as determined from molecular dynamics computer simulations. 
Furthermore we calculate the wave vector dependence of the Debye- Waller 
factor for a realistic model of silica and compare these results with the ones 
obtained from a simulation of this model. We find that if the contribution 
of the three point correlation function to the vertices of the memory kernel 
are taken into account, the agreement between theory and simulation is very 
good. Hence we conclude that mode coupling theory is able to give a correct 
quantitative description of the caging phenomena in fragile as well as strong 
glass-forming liquids. 
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I. INTRODUCTION 

In the last decade our understanding of the dynamics of supercooled liquids has made 
significant progress In particular it has been shown that for fragile glass- formers the 
bend one observes if one plots the logarithm of the viscosity as a function of the inverse 
temperature can be explained very well by means of the so-called mode-coupling theory of 
the glass transition (MCT) 0. In the vicinity of this bend the dynamics of the system 
changes qualitatively in that the particles start to experience strong caging effects, i.e. they 
are temporarily trapped by the particles that surround them. MCT gives a self-consistent 
description of the dynamics of the particles inside this cage as well as how the particles leave 
this cage, i.e. of the structural relaxation of the supercooled liquid. In the past the predic- 
tions of this theory have been checked in many experiments as well as computer simulations 
and it was found that MCT is indeed able to give a qualitatively correct description of the 
relaxation dynamics p|. 

However, in principle the theory is supposed to give not only a qualitatively description 
of the relaxation dynamics of supercooled liquids, but also a quantitative one, if the static 
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properties of the system are known with sufficiently high precision. This attractive feature 
originates directly from the way the theory is (or can be) derived, namely the Mori-Zwanzig 
formalism in which one obtains equations of motion for slow variables which involve their 
static values. Thus, once these static values are known one can, in principle, determine their 
time dependence. In particular it is possible to calculate from the knowledge of the static 
structure factor the time dependence of the coherent and incoherent intermediate scattering 
functions, F{q,t) and Fs{q,t), respectively, where q is the wave vector. Unfortunately, in 
the past these type of calculations have been done only for very few systems, since one the 
one hand they are quite involved and on the other hand they require as input structural 
data with very high quality (better than 1%) 1^-0. In the present paper we expand these 
type of calculations in two directions. On the one hand we solve for a simple glass-forming 
system, a binary Lennard- Jones mixture (BMLJ), the full time and wave vector dependence 
of the MOT equations, including a realistic short time dynamics, and compare them with 
results from computer simulations of the same system. On the other hand we calculate 
the g— dependence of the nonergodicity parameters (NEP) for silica (Si02), a glass-former 
whose structure is given by an open tetrahedral network and who is the prototype of a 
strong glass-former, and compare also these results with the ones from simulations of the 
same system. 

II. THEORY 

In this section we summarize the MCT equations that are needed to calculate the quan- 
tities discussed in section In order to keep this presentation as simple as possible we 
will discuss only the equation for the case of a one-component system, although in reality 
we have used the equations for a two-component system, since the BMLJ as well as Si02 
belong to this class. The full binary equations can be found in Refs. [3, P^ , P^] • 

The intermediate scattering function can be defined by F{q, t) = (5p(q, 0)(5p(q, t)) where 
6p are the density fluctuations. F{q,t) obeys the exact equation of motion 

F(g, t) + n\q)F{q, t) + f M{q, T)F{q, t - T)dT = (1) 

where the frequency Vl^ is given by Vt^ = q^ksT / {mS{q)) . Here m and S{q) are the mass of 
the particles and the static structure factor respectively. The function M(g, t) in Eq. (|ID is 
the so-called memory function and it is useful to write it as follows: 

M(g,t) = M^^«(g,t) + {m*^^^[F(A;, t)](g) - M*^^^[F^(fc, t)](g)} . (2) 

Here M^^^{q, t) is that part of the memory function which is responsible for the dynam- 
ics of the system at very short times, i.e. after the particles have left the ballistic regime. 
The functional M^^^[F(fc,t)] is the usual memory kernel of MCT which depends on the 
static structure factor as well as on the three point correlation function C3(q, k) 0. Also 
it contains a short time part. But since we want to describe this time regime by means of 
M^^'^{q,t), we have to subtract out this part from M^^*^^[F(/c, t)]. This is done in the last 
term of Eq. (^, where F^ is a function which decays rapidly to zero, but has the correct 
behavior at short times [p!0|] . 
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For the regular memory function M^^'^{q,t) one can make different type of Ansatzes. 
One which seems to work well at high temperatures is given by [|l^ 

M^^«(g,t) = a(g)/cosh(/3(g)t). (3) 

Here a{q) and (3{q) are parameters which can be calculated via sum rules from the static 
structure factor and other static quantities which can be measured in a computer simula- 
tion ||T0|jr3[| . Hence they are not adjustable fit parameters. 

Eqs. are a self-contained set of equations of motion from which one thus can 

calculate the full time and wave vector dependence of F{q, t), and similar equations exist for 
the incoherent intermediate scattering function Fs{q,t). We have solved these equations by 
discretizing g-space into 100 points that covered the g— range up to 3-4 times the location 
of the main peak in S{q). 



III. MODELS AND DETAILS OF THE SIMULATIONS 

The first model investigated is a 80:20 mixture of Lennard- Jones particles. In the fol- 
lowing we will call the majority and minority species A and B particles, respectively. Both 
of them have the same mass m and they interact via a potential Vajs = ^€.afi[{o'ap/'rY'^ — 
{o'ai3/r)^], a,P & {A, B}. The parameters and cTq,/? are given by caa = 1-0, ctaa = 1-0, 
eAB = 1-5, (Tab = 0.8, cbb = 0.5, and ctbb = 0.88. This potential is truncated and shifted at 
a distance aa/3- In the following we will use cr aa and eAA as the unit of length and energy, 
respectively (setting the Boltzmann constant = 1.0). Time will be measured in units 



of yma\j^/ ASeAA- In the past the structural and dynamical properties of this system have 
been studied in great detail [|14|,|15|. More detail on this can be found in Ref. [|16 . 



For the present work we only needed to determine the three point correlation function cs 
since the time and temperature dependence of F{q, t) and Fs{q, t), as well as the one of S{q), 
can be found in the mentioned literature. For this we simulated a system of 800 A particles 
and 200 B particles in a box with volume (9.4)'^. The total time of this simulation was 
about 10^ time steps from which we obtained roughly 12,000 independent configurations. 
This large number was necessary to determine C3 with sufficient precision. Due to this large 
computational effort we did this calculation only for one temperature, T = 1.0. Thus in the 
following we will assume that the temperature dependence of C3 is weak. 

The second model we study is amorphous silica, Si02. For this we use the potential 
proposed by van Beest et al. which has the functional form |T7| 



/.^^(^) = M^ + ^^^exp(-i?,^r)-^ «,/9e[Si,0]. (4) 



The values of the constants qa,qf3, Aajj, Baf3, and C^/? can be found in Ref. |T^. The 
potential has been truncated and shifted at 5.5 A. In the past it has been shown that this 
potential is able to give a reliable description of silica in its molten phase as well as in the 



glass (see PB|-[2D[] and references therein). For the present calculations to determine C3 we 
used 600 ions in a box with volume (20. 4A)^. The total length of the simulation was 2 x 10^ 
time steps, from which we obtained at 4000 K around 2000 independent configurations. 



3 



IV. RESULTS 



We start by considering first the dynamics of tlie BMLJ at intermediate and liigfi tem- 
peratures. In tliis T— range it can be expected tliat tlie effect of tlie memory kernel of MCT 
is not relevant and thus we will set it to zero. In Fig. |I| we show the time dependence of 
Faa{(1, t) for various temperatures. The wave vector is 7.25, the location of the main peak in 
Saa{<i)- The dashed lines with symbols are the result of the simulation whereas the full lines 
are the prediction of the theory. As can be seen, the theory works very well at high tem- 
peratures but starts to break down at intermediate temperatures in that it underestimated 
the correlation function at intermediate times. Thus we see that even at the intermediate 
temperature T = 1.0, which is more than twice the MCT temperature Tc = 0.435, cage 
effects become important. 

In order to see whether the memory kernel M'^^'^ is able to take into account these 
effects we have solved Eqs. by taking now into account also this contribution to the 

memory function M(g, t). In doing this we had to face a problem which we had encountered 
already some time ago 0, namely that MCT is not able to predict reliably the value of the 



critical temperature Tc. For the BMLJ the simulations show that Tc ~ 0.435 |T^, whereas 
the theory predicts a value around 0.92 This means that the theory is not able to predict 
correctly the absolute value of the time scale for the a— relaxation, although it is able to 
predict the shape of the correlation functions (see below). Therefore we had to use one 
adjustable parameter, a temperature which we will denote by Tj, which is the temperature 
at which the vertices in the MCT-functional M^^'^ are evaluated. The value of Tj was 
adjusted such that the time scale for Taa(q', ^) for q = 7.25 was reproduced correctly. 

In Fig. ^ we show the time dependence of the coherent as well as the incoherent inter- 
mediate scattering function for q = 7.25 and q = 9.98, the location of the first peak and the 
first minimum in S'aa(q'); respectively. The temperature is 2.0, i.e. a value for which we find 
that the regular memory kernel is no longer able to give a good description of the relaxation 
dynamics (see Fig. From this figure we see that in general the agreement between the 
simulation and the theory is very good in that the shape of the curves as well as their po- 
sition is correctly predicted. (The discrepancy found for Faa{(1 = 9.98, t), where the theory 
predicts a pronounced shoulder at around t = 2 whereas the simulation shows only a weak 
shoulder in that time regime, is probably related to the fact that the description of the short 
time dynamics is not yet completely adequate ||10[)- Thus we conclude that MCT is indeed 
able to give a correct description of the relaxation dynamic of the system at intermediate 
temperatures, i.e. at temperatures where the cage effect starts to become noticeable. 

We now check whether this conclusion is also correct if the temperature is so low that 
the cage effect becomes very important. For this we have solved the MCT equations for 
T = 0.466, i.e. a temperature for which the relaxation dynamics is about lO'^ times slower 
than the one at high T. The time dependence of -Faa(q', t) and of Fl{q, t) as predicted from 
the theory is shown in Fig. |[ Also included are the results from the computer simulations 
from Ref. As in the case of intermediate temperatures. Fig. |^, we find that also for 

this T the agreement between theory and simulation is very good. The main discrepancy is 
again seen for T'aa(q', t) aX q = 9.98, and the reason for it is the same as the one given above. 
All in all we thus conclude that for this system the theory is indeed able to predict the full 
time and wave vector dependence of the coherent and incoherent scattering function. 
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The temperature dependence of the relaxation time of the BML J system shows significant 



deviations from an Arrhenius law [14|. As mentioned above, these deviations are believed 



to be related to a change in the transport mechanism of the particles which show a hopping 
type of motion at low temperatures whereas at high T they show a more collective/flow- 
like behavior. It is of interest that recently it has been suggested that even silica shows 
such a crossover in the transport mechanism, although this crossover occurs at relatively 
high temperatures (around 3300K) |TB[. Therefore one might ask whether MCT is able to 
give a reliable description of the relaxation dynamics of this important glass-forming system 
also. Note that from a structural point of view the BMLJ system and amorphous silica are 
very different, since the former one resembles to the random close packing of hard spheres 
whereas the latter is given by an open network structure similar to the continuous random 
tetrahedral network proposed long time ago by Zachariasen Since, as pointed out in 
the Introduction, MCT uses only structural information to predict the dynamics, it is of 
great interest to see whether the theory is also able to give a correct quantitative description 
if the structure is very different from the one of closed packed hard spheres. 

To check this we have calculated for amorphous silica the wave vector dependence of 
the nonergodicity parameters, i.e. the height of the plateau in the intermediate scattering 
function at intermediate times (see, e.g.. Fig. ^. Before we discuss the results, we have to 
mention a technical point which makes the calculation of the NEP for the case of silica much 
harder than for the case of the BMLJ. In Sec. H we mentioned that the memory kernel of 
the MCT contains only static quantities, namely the static structure factor S{q) and the 
three point correlation function C3(q, k). From a simulation it is quite easy to determine 
S{q) with high precision. For the function C3 this is, however, not the case, since due to the 
two vectorial arguments the statistics for this quantity is very bad. Therefore we had to 
make very long simulations in order to determine C3 with sufficient accuracy. More details 
on this can be found in Ref. [|ll|]. 

In the following we will discuss the results for the NEP for the BMLJ as well as for 
the case of silica. Since all of the results presented in Figs. |l|-0 were obtained with the 
approximation that C3 = 0, one has of course to check whether or not they do not change if 
this assumption is not made. We mention, however, already here, that some time ago Barrat 
et al. showed that this approximation is very good for the case of a soft-sphere system, i.e. 



a system which is relatively similar to the BMLJ considered here p2[. Whether this result 



holds also for the case of a system with an open network structure has, however, so far not 
been investigated. We also mention that in order to calculate the NEPs it is not necessary to 
introduce any fit parameter of any kind. The only input to the data are the partial structure 
factors [^. Thus for this type of calculation the above discussed problem with the Tf does 
not exist. 

In Fig. ^ we show the wave vector dependence of the NEP for the coherent functions. 
(Note that since this is a binary system, there are three of them. Furthermore we mention 
that for reasons of convenience we show the NEP multiplied by the corresponding partial 
structure factors.) In each panel we show three curves: The circles are the result from the 
simulation published in Ref. The dashed and full line is the theoretical result for the 
cases that C3 is set equal to zero and C3 7^ 0, respectively. First of all this figure shows that 
the theory is able to reproduce with excellent accuracy the data of the simulation without 
any adjustable parameter. Furthermore we recognize from Fig. 11 also that the theoretical 
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prediction hardly depends on whether or not C3 is taken into account, in agreement with the 
finding of Barrat et al. ^ . 



For sihca the situation is quite different as can be inferred from Fig. ^ where we show the 
wave vector dependence of the NEP for this system. We see that in this case the theoretical 
prediction for C3 = differs strongly from the one if this function is taken into account. 
Thus we find that for the case of a network structure the contribution of the three point 
correlation function to the memory function is very important. It is remarkable that if the 
contributions of C3 are taken into account, the theoretical prediction agrees very well with 
the result of the simulation, which were presented in Ref. [rU|. Thus we conclude that the 
theory is also able to give a quantitative correct prediction for this type of glass-former. 



V. SUMMARY 

The goal of this work was to check to what extend the mode-coupling theory of the 
glass transition is able to give a correct quantitative description of the dynamics of glass- 
forming liquids. This was done for two very different types of systems: A binary mixture 
of Lennard- Jones particles, whose structure is similar to the one of a close packing of hard 
spheres and whose temperature dependence of relaxation times makes it a glass-forming 
liquid of intermediate fragility. On the other hand we have studied silica, who has a open 
network structure and which, in the temperature region where experiments are feasible, is 
considered to be the prototype of a strong glass-former. 

For the BMLJ system we have solved the MCT equations to obtain the full time and 
wave vector dependence of the relaxation dynamics. In particular we have also included a 
realistic Ansatz for the dynamics at short times so that the theoretical curves should be able 
to give also a good description in this time regime. By comparing the theoretical curves for 
F{q,t) and Fs{q,t) with the ones obtained from computer simulations of the same system, 
we find that cage effects become noticeable already at relatively high temperatures. The 
theory is able to give a very reliable quantitative description of the relaxation dynamics for 
all temperatures considered. The only discrepancy found is probably related to the fact that 
our understanding of the dynamics at short times is still incomplete. 

For the case of silica we have calculated the g— dependence of the nonergodicity param- 
eters. We have found that for this system it is important to include in the evaluation of the 
memory function also those contributions that stem from three point correlation functions. 
Most probably this finding is related to the open network structure of the system. We find 
that once C3 is taken into account the prediction of the theory for the NEP are in very good 
agreement with the results from computer simulations. 

In summary we have shown that MCT is able to give a very good quantitative descrip- 
tion of the relaxation dynamics of fragile as well as strong liquids, at least at high and 
intermediate temperatures. Thus we conclude that this theory is able to rationalize at least 
the first few decades of the slowing down of a very large class of glass-forming liquids on a 
quantitative level. 
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FIGURES 




FIG. 1. Time dependence of the coherent intermediate scattering function for the A particles in 
the BMLJ system for different temperatures as obtained from the simulation (dashed lines with 
symbols) and as predicted from the theory if the MCT kernel is not taken into account (solid lines). 
For clarity some of the curves have been shifted vertically. 
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FIG. 2. Time dependence of the coherent and incoherent intermediate scattering function of the 
BMLJ system for two wave vectors at T = 2.0. The dashed line with the symbols arc the results 
from the simulation and the solid lines are the prediction of the theory in which the MCT kernel 
has been taken into account. 
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FIG. 3. Same as figure 3 but now for T = 0.466. 
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FIG. 4. Wave vector dependence of rh#Tionergodicity parameters for the BMLJ system. The 
circles are the result of the simulation, the dashed line is the theoretical prediction if the three 
point correlation function is set to zero, and the full line is the theoretical prediction if this function 
is taken into account. 
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